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BALLISTIC  RESEARCH  LABORATORIES 

MEMORANDUM  REPORT  NO.  6l? 


WBFritz/graw 

Aberdeen  Proving  Ground,  Md. 
August  1952 


A  SURVEY  OF  ENIAC  OPERATIONS  AND  PROBLEMS: 

19U6  -  1952 


ABSTRACT 


This  report  summarizes  the  entire  operating  history  of  the  Eniac 
since  Fdnruary  I9U6,  when  the  machine  was  demonstrated  to  the  public 
as  the  first  large  scale  general  purpose  electronic  digital  computing 
machine.  It  brings  up  to  date  much  of  the  material  presented  in  the 
unpublished  paper  entitled  wThe  Eniac  -  A  Five  Year  Operating  Survey” 
presented  at  the  Detroit  meeting  of  the  Association  for  Computing 
Machinery  on  27  March  1951* 

The  primary  interest  of  the  paper  is  the  brief  descriptions  given 
of  the  specific  unclassified  problems  completed  during  this  period. 

It  is  hoped  that  individuals  not  yet  in  immediate  contact  with  high 
speed  computers  will  gain  some  insight  into  the  type  of  problems  which 
machines  make  accessible  to  computation.  Pertinent  reference  material 
is  given  to  a  variety  of  unclassified  published  material  related  to 
these  couplet ed  programs. 
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Intentionally  left  blank. 


I.  HISTORICAL  INTRODUCTION 


The  Eniac  -  the  Electronic  Numerical  Integrator  and  Computer  - 
was  the  first  high  speed,  electronic,  general  purpose,  digital  com¬ 
puting  machine*  It  was  designed  and  constructed  at  the  Moore  School  of 
Electrical  Engineering  of  the  University  of  Pennsylvania  (l^ljlj.  -  19^6). 

The  primary  reason  for  the  construction  of  the  Eniac  was  to  speed  up 
the  large  quantity  of  numerical  computations  required  for  the  Army's 
firing  tables  and  bombing  tables*  Since  the  Eniac  was  designed  to 
handle  this  particular  type  of  problem,  there  was  little  attempt  to 
simplify  the  steps  necessary  to  prepare  new  problems  or  to  change  from 
one  problem  to  another*  The  original  method  of  Eniac  control  known  as 
"direct  programming"  required  the  activation  of  a  large  number  of  cir¬ 
cuits  by  the  tedious  procedure  of  setting  switches  and  connecting  plug 
wires  in  some  forty  panels  in  addition  to  setting  tabulated  data  on 
the  switches  of  the  "portable"  function  tables,  whenever  a  problem  was 
to  be  put  on  the  machine*  The  size  of  the  problem  that  could  be  solved 
was  limited  by  the  number  of  controls  on  each  of  the  machine's  units* 

The  number  of  problems  that  could  be  handled  within  a  given  period  of 
time  was  limited  by  the  time  it  took  to  prepare  the  program  of  control 
settings  necessary  to  solve  the  problem. 

Experience  with  the  first  problems  solved  on  the  Eniac  at  Philadelphia 
and  later  at  Aberdeen,  showed  clearly  that  the  cumbersome  method  of  direct 
programming  would  severely  limit  the  Eniac 's  use  as  a  general  purpose 
computer.  Operating  percentages  were  not  kept  during  this  early  period, 
but  various  estimates  of  machine  productivity  have  been  made  ranging 
from  to  as  much  as  2$%a  However,  during  this  period,  the  Eniac  was 

only  in  operation  for  a  normal  8  hour  shift  for  $  days  each  week*  All 
voltages  except  filament  voltages  were  turned  off  during  the  non-operative 
period.  To  take  advantage  of  a  productive  machine,  however,  the  oper¬ 
ator  occasionally  overstayed  the  normal  I63O  departure  time.  During 
this  period  of  the  machine's  operation,  from  February  19U6  to  March  19^8, 
eleven  problems  were  couplet ed. 

Early  in  19^8,  R.  F.  Clippinger  and  some  of  his  associates,  in  the 
course  of  coding  the  solution  of  a  sinple’  system  of  five  simultaneous 
hyperbolic  partial  differential  equations,  were  forced  to  adopt  a  dif¬ 
ferent  method  of  using  the  Eniac  in  order  to  fit  their  problem  on  the 
machine.  This  method  consisted  of  a  code  for  parallel  operation  using 
numerical  information  from  its  function  tables  for  logical  control.  The 
experience  with  this  method  (first  discussed  in  reference  1),  led 
J.  von  Neumann  to  suggest  the  use  of  a  serial  code  for  control  of  the 
Eniac.  Such  a  code  was  devised  and  employed  with  the  Eniac  beginning 
in  March  191*8.  Operation  of  the  Eniac  with  this  code  was  several  times 
slower  than  either  the  original  method  of  direct  programming  or  the  code 
for  parallel  operation.  However,  the  resulting  siuplification  of  coding 
techniques  and  other  advantages  far  outweighed  this  disadvantage. 


II.  BRIEF  DESCRIPTION  OF  THE  ENIAC 


The  Eniac  is  at  present  controlled  by  an  improved  version  of  the 
original  one  address  type  serial  code,  which  incorporates  a  unit  called 
the  Converter  as  a  basic  part  of  its  operation,  hence  the  name,  the 
Eniac  Converter  Code.  The  continual  improvement  of  this  code  has  been 
brought  about  primarily  by  changes  in  coding  techniques  and  engineering 
additions  to  the  original  machine.  The  Eniac,  however,  is  still  very 
limited  in  high  speed  memory  and  somewhat  awkward  in  comparison  to  the 
present  day  concept  of  the  capabilities  of  a  large  scale  digital  com¬ 
puter.  The  Eniac  has  only  fifteen  high  speed,  storage  registers  in 
addition  to  the  five  required  primarily  as  part  of  the  arithmetic  unit. 
Each  register  or  accumulator  is  able  to  store  a  signed  ten  digit  decimal 
number. 

Of  the  2i;00  individual  2  digit  instructions  of  a  single  program, 
l821j.  must  be  set  by  hand  on  the  36U8  switches  of  the  control  panels. 

The  remaining  576  instructions  are  wired  on  IBM  plugboards  which  are 
prepared  prior  to  the  actual  machine  set  up  period.  The  average  time 
required  to  place  a  problem  on  the  Eniac  is  about  1.5  hours. 

The  prime  advantage  of  the  Eniac  is  its  speed.  The  usual  single 
instruction  arithmetic  operations  are  performed  at  the  following  ap¬ 
proximate  rates:  1)  Addition  -  l600/sec,  2)  Subtraction  -  1000/sec, 

3)  Multiplication  -  290/sec,  i;)  Division  and  Square  Rooting  -  i|0/sec. 

Of  course  in  one  address  type  codes,  a  single  arithmetic  operation  may 
require  more  than  one  instruction  depending  on  the  availability  of  the 
quantities.  If  we  assume  the  Eniac  were  to  perform  only  isolated  oper¬ 
ations  unrelated  to  the  rest  of  the  computation,  the  following  operation 
rates  would  be  more  correct:  l)  Addition  -  550/sec,  2)  Subtraction  - 

U50/sec,  3)  Multiplication  -  192/sec,  and  U)  Division  and  Square 
Rooting  -  32/sec.  Actual  production  figures  for  efficient  coding, 
however,  show  rates  more  closely  approximating  the  first  set  of  figures. 

The  term  approximate  rate  has  been  used  since  a  variable  oscil¬ 
lator  controls  the  frequency  of  the  pulse  rate,  hence  the  speed  of  the 
operations.  Prior  to  the  adoption  of  the  code,  the  Eniac  was  always 
operated  at  a  frequency  of  100  KC.  During  the  early  period  of  the 
code  (19l;8-19H9),  the  operating  frequency  was  often  as  low  as  60-70  KC. 

At  the  present  time,  the  average  rate  is  85  KC. 

Another  advantage  of  the  Eniac  is  the  efficient  input-output  de¬ 
vices  -  i.e.  an  IBM  reader  and  an  IBM  punch.  The  8  signed  10  digit 
numbers  on  an  IBM  card  can  be  transferred  to  the  high  speed  memory  in 
about  0.3  seconds  or  essentially  at  the  rate  of  26  words  per  second. 

A  similar  group  of  80  decimal  digits  from  the  high  speed  memory  can 
be  punched  on  an  IBM  card  in  0.2  to  0.6  seconds  depending  on  the  fre¬ 
quency  of  the  punching.  This  means  a  rate  of  13  to  HO  words  per  second, 
with  the  usual  rate  being  1±0  per  second.  Tabulated  data  or  constants 
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may  be  brought  into  the  electronic  memory  from  the  hand  set  control 
panels  at  the  rate  of  715  quantities  per  second.  Certain  problems  re¬ 
quiring  a  large  memory  have  been  run  efficiently  on  the  Eniac  by  punch¬ 
ing  intermediate  results  on  IM  cards  which  in  turn  were  read  back  into 
the  machine  for  further  computation.  A  complete  description  of  the  9k 
arithmetic,  storage,  and  control  instructions  of  the  Eniac  Code  is 
given  in  Reference  2. 


III.  OPERATING  EXPERIENCE 


Experience  over  the  past  several  years  has  demonstrated  that  the 
code  has  improved  the  productive  operation  of  the  Eniac  in  the  following 
ways :  It 

1*  made  easier  the  preparation  of  problems  for  the  machine, 

2.  permitted  an  increase  in  the  size  of  problem  that  can  be 
handled  by  the  machine, 

3.  reduced  the  time  required  for  changing  problems  from  several 
days  to  several  hours,  and 

i;.  simplified  testing  procedures. 

At  the  time  the  Converter  Code  was  added  to  the  Eniac,  J.  V.  Holberton 
initiated  a  system  of  record  keeping  which  accounts  for  every  minute  of 
the  time  the  machine  is  manned.  This  log  of  the  machine’s  operation  has 
been  continuously  kept  since  April  19^8.  The  quarterly  records  have  been 
summarized  in  Figure  1. 

The  increase  in  trouble-shooting  time  since  1  April  1951  is  par¬ 
tially  explained  by  the  transfer  of  several  key  service  engineers,  a 
general  emphasis  on  the  part  of  the  Computer  Research  Branch  to  complete 
the  Edvac,  the  failure  to  obtain  certain  required  replacement  parts, 
and  a  general  upswing  in  the  number  of  problems  handled  by  the  Eniac. 

Some  of  the  obvious  requirements  for  good  machine  operation  are 
coupe  tent  operating  and  engineering  personnel,  efficient  testing  pro¬ 
cedures  and  adequate  programming.  Another  factor  affecting  operating 
percentages  is  the  degree  of  human  intervention  which  is  used  in  solving 
the  problem.  Frequently  no  human  intervention  is  required.  However, 
one  type  of  problem  may  require  that  the  mathematician  direct  the  course 
of  the  computation  by  determining  each  new  step  on  the  basis  of  the  re¬ 
sults  of  the  previous  steps.  If  possible  this  type  of  machine  usage 
is  avoided  but  occasionally  it  is  not  possible  to  state  numerically, 
ahead  of  time,  precise  criteria  for  various  courses  of  action.  Data 
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ENIAC  CONVERTER  CODE¬ 
MACHINE  EFFICIENCIES  PER  QUARTER 


00 


% 


1.  O - O  CORRECTLY  OPERATING  ON  THE  SOLUTION  OF  REGULAR  PROBLEMS. 

2.  □ - □  LOCATING  AND  CORRECTING  MACHINE  TROUBLE 

IN  THE  ENIAC,  NON  DUPLICATION  TIME,  AND  DOWN  TIME  ON  SPECIAL 
PREVENTIVE  MAINTENANCE. 

3  A— - A  PLACING  NEW  PROBLEMS  ON  THE  ENIAC,  CHECKING 

PROGRAMMING,  DATA  ANALYSIS,  AND  DOWN  TIME  DUE  TO  HUMAN 
OPERATING  ERROR. 


FIGURE  I 


analysis  time  in  a  problem  of  this  type  will  be  large.  Another  type 
of  problem  that  of  necessity  results  in  reduced  operating  efficiency 
is  one  in  which  the  control  panels  have  to  be  changed  during  the  course 
of  the  problem.  Here  the  memory  is  insufficient  to  handle  the  problem 
efficiently;  but  the  Eniac,  at  the  time,  was  the  only  means  available 
to  obtain  the  desired  computations  in  a  reasonable  period  of  time. 

The  figures  indicated  in  Figure  1  are  percentages.  The  actual 
records  show  an  average  of  over  l£00  machine  hours  per  quarter  (approxi¬ 
mately  116  hours  per  week).  The  normal  work  week  has  been  0800  Monday 
through  2I4.OO  Friday. 

Putting  this  productivity  another  way,  the  machine  has  operated 
correctly  on  the  solution  of  problems  of  the  time  for  the  period 
1  April  I9I4.8  through  30  June  19  £2. 

It  is  interesting  to  note  that  a  staff  of  l£  -  2£  mathematicians, 
technicians,  and  engineers  have  been  utilized  to  perform  programming, 
service  engineering,  and  to  operate  the  aiiac. 

Before  giving  some  information  on  the  number  and  type  of  problems 
handled  by  the  machine,  it  might  be  well  to  define  just  what  we  mean  by 
a  "problem".  A  "problem"  is  defined  as  a  machine  program  or  group 
of  machine  programs  related  to  a  specific  subject  under  investigation 
or  a  particular  mission  to  be  accomplished.  Ex  apples  of  single  problems 
are  the  following: 

1.  The  Doppler  determination  of  piissile  position.  Here  the 
single  "problem"  can  be  adapted  to  any  required  number  of 
individual  firings,  hence  the  problem  is  repeated  periodically 
with  different  input  data.  Various  improvements  have  been 
made  and  the  program  of  instructions  increased  to  make  mor$ 
autopiatic  the  determination  of  the  required  solutions. 

2.  Determination  of  the  mole  fractions  of  the  gaseous  con¬ 
stituents  of  the  equilibrium  composition  of  four  component 
systems  of  hydrogen,  oxygen,  nitrogen,  and  carbon  at  tem¬ 
peratures  of  1000°K  to  £000®K  and  pressures  of  0.1  atmos¬ 
pheres  to  those  of  interior  ballistics.  In  this  instance 

a  number  of  sets  of  equations  and  a  large  range  of  parameters 
are  involved.  Several  machine  months  have  been  expended  on 
the  problem  in  a  number  of  short  machine  periods. 

A  large  portion  of  the  problems  completed  on  the  Eniac  are  dis¬ 
cussed  in  the  next  section. 
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17.  GUIDE  TO  COMPLETED  PROBLEMS 


Since  the  Eniac  was  completed  in  February  19U6,  nearly  100  problems 
have  been  completed.  It  will  be  impossible  to  go  into  much  detail  de¬ 
scribing  these  problems,  but  it  is  felt  that  a  brief  guide  to  the  work 
that  has  been  accomplished  would  be  of  general  interest. 

Some  70%  of  these  problems  have  originated  within  the  Ballistic 
Research  Laboratories,  while  the  remainder  were  submitted  by  other 
government  agencies,  or  either  private  industrial  organizations  or 
universities  with  government  contracts.  Over  half  of  the  problems 
have  had  to  do  with  the  numerical  int egration  of  systems  of  ordinary 
differential  equations.  The  following  integration  methods  are' typical 
and  have  been  used  to  solve  one  or  more  problems:  Heun,  Euler,  Runge- 
Kutta,  Iterated  Simpson,  and  Kutta-Sinpson.  One  fifth  of  the  problems 
have  been  primarily  solutions  of  algebraic  systems.  Of  equal  frequency 
has  been  the  solution  of  partial  differential  equations  including  sys¬ 
tems  of  hyperbolic  equations  and  single  equations  of  the  parabolic 
type.  The  remaining  problems  have  been  the  solution  of  integral 
equations  or  number  theoretic  computations. 

The  next  few  pages  survey  a  large  group  of  problems  completed  on 
the  Eniac,  the  problems  being  listed  in  connection  with  the  group 
sponsoring  the  computation.  It  should  be  pointed  out  that  security 
regulations  prohibit  a  complete  listing  of  all  the  work  performed. 

If  specific  unclassified  Eniac  problems  with  which  the  interested 
reader  is  familiar  are  not  mentioned,  it  is  suggested  that  he  contact 
the  author  in  order  that  any  future  "guides1'  may  correct  these  over¬ 
sights. 

The  following  chart  summarizes  the  work  discussed  in  the  remainder 
of  the  report. 


Sponsor 

No.  of  Problems 

%  Machine  Time 

I 

Computing  Laboratory 

35 

II. 

Exterior  Ballistic  Lab. 

8 

12 

Ill 

Interior  Ballistic  Lab. 

3 

2 

17 

Ballistic  Measurements  Lab. 

1 

1 

V 

Surveillance  Laboratory 

1 

1 

71 

Outside  Groups 

20 

h9 
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I*  Computing  Laboratory 

The  largest  amount  of  computing  for  any  single  group  was  naturally 
enough  performed  for  the  Computing  Laboratory  of  BRL.  These  confutations 
are  discussed  in  three  categories: 

A«  Ballistic  Tables 

B*  Measurements  Analysis 

C.  Various  Mathematical  Investigations 


A*  Ballistic  Tables 

In  recent  years  about  of  the  Eniac's  computing  time  has  been 
spent  on  the  development  of  Ballistic  Tables  involving  confutations  on 
some  lil  problems.  The  determination  of  the  numerical  data  for  Firing 
Tables  and  Bombing  Tables  is  a  more  extensive  conp utational  problem  than 
the  mere  confutation  of  the  trajectories.  The  following  mathematical 
problems  have  been  programmed  in  order  to  determine  the  final  form  of 
those  tables: 

1.  Solution  of  the  non-linear  differential  equations 
describing  the  path  of  a  trajectory. 

2.  Enpirical  determination  of  parameters  by  least  square 
techniques  to  adapt  theory  to  practice,  and 

3.  Extensive  high  order  interpolation  in  order  to  present 
the  results  in  usable  form. 

A  brief  survey  of  some  of  the  Eniac's  work  on  ballistic  tables 
is  given  in  "Proceedings  of  the  Pittsburgh  Meeting  of  the  ACM  2-3  May 
1952."  "Firing  Table  Computations  on  the  Eniac"  -  H.  L.  Reed,  Jr. 

B.  Measurements  Analysis 

The  Eniac  has  been  utilized  in  the  reduction  of  data  from  high 
speed  cameras  in  order  to  obtain  certain  information  regarding  the 
flight  of  guided  and  unguided  missiles  and  bombs. 

\ 

The  reduction  problem  to  determine  the  azimuth,  elevation, 
yaw,  and  orientation  of  yaw  of  the  missile  axis  of  bombs  and  missiles 
is  described  in  the  following  reports: 

1.  BRL  Ilk  -  "Attitude  and  Yaw  Reductions-  of  Projectiles  in 

Free  Flight"  -  G.  R.  Trimble,  Jr.  -  19J>1. 

2.  BRL  812  -  "A  Least  Square  Attitude  Solution"  - 

G.  R.  Trimble,  Jr.  -  1952* 
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The  determination  of  space  coordinates,  and  velocity  and  ac¬ 
celeration  components  of  ground  to  ground  missiles  as  well  as  bomb 
releases  are  obtained  from  high  speed  camera  data.  Observations  from 
two  or  three  camera  stations  provide  azimuth  and  elevation  angles  of 
the  line  of  sight,  A  least  square  procedure  is  used  to  minimize  the 
distance  between  the  lines  of  sight,  from  which  the  space  coordinates 
are  readily  obtained. 

3,  BRIM  $l\3  -  "The  Preparation  of  a  Theodolite  Problem  for 
Large  Scale  Computing  Machines"  -  B.  Dimsdale 
and  M.  A*  Powers  -  1951* 

h»  BRL  778  -  "A  Recursion  Relation  for  Computing  Least- 

Square  Polynomials  over  Moving  Arc"  - 
G.  R.  Trimble,  Jr,  -  1951* 

C,  Various  Mathematical  Investigations: 

The  following  group  of  problems  concern  a  variety  of  mathe¬ 
matical  investigations  as  well  as  number  theoretic  computations. 

The  work  of  several  individuals  closely  associated  with  the 
Computing  Laboratory  is  included  here, 

1,  Theory  of  Rounding  Errors 

An  investigation  of  error  accumulation  was  performed 
and  reported  on  by  Hans  Rademacher  as  Lecture  19  in  the 
series  of  lectures  given  at  the  Moore  School  of  the 
University  of  Pennsylvania  -  8  July  -  31  August  19^6, 

Hans  Rademacher' s  paper  was  entitled  "On  the  Accumulation 
of  Errors  in  Numerical  Integration  on  the  Eniac"  and  was 
published  in  Volume  II  of  the  four  volume  set  on  the 
lecture  series,  pp  19-1  —  19-17*  The  equations  integrated 
were  y*  *  p  and  p  *  -y  thus  generating  sines  and  cosines. 

2,  Mil  Tables  of  Sines  and  Cosines 

Mil  tables  of  the  sine  and  cosine  for  0(1)  1600  to 
10  decimal  places  using  the  formula  cos(n  +  1)  0  ■ 

2  cos  0  cos  n  8  -  cos  (n-1)  0,  The  table  was  repro¬ 
duced  locally  for  use  in  ballistic  work. 

3*  Investigations  in  Prime  Numbers 

The  confutations  completed  during  several  holiday 
week  ends  by  D.  H,  Lehmer  of  the  University  of  California 
have  been  reported  in  three  papers. 
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a.  Math.  Tables  and  Other  Aids  to  Computation  v.  2, 

p.  313,  '1*7. 

b.  Amer.  Math.  Soe.  Bulletin  v.  53,  19l±7,  pp  I6I1-I67,  "On  the 

Factors  of  2n  +.  1"  by  D.  H.  Lehmer. 

c.  The  American  Mathematical  Monthly  v.  56,  No.  5, 

pp  300-309,  Kay  1 9h9 ,  "On  the  Converse  of 
Fermat’s  Theorem  H"  -  D.  H.  Lehmer. 

it.  Fermat’s  Quotient 

Determination  of  a  Table  of  Fermat’s  Quotients  obtained 
during  machine  check  periods.  MTAC  v.  5,  Apr.  1951, 

pp  81|-85. 

5.  Calculation  of  u,  e,  and  1/e 

u  was  computed  to  2035  places  and  e  to  2556  places 
in  an  attempt  to  gain  insight  into  the  character  of 
these  numbers  by  obtaining  the  distribution  of  a  large 
number  of  their  digits  as  well  as  an  illustration  of 
machine  methods  for  handling  a  large  number  of  decimal 
places.  1/e  was  computed  as  a  check  on  e.  e  *  1/e  was 
obtained  on  an  IBM  tabulator. 

a.  MTAC  17  No.  29  -  Jan.  1950,  "An  Eniac  Determination 

of  n  and  e  to  More  than  2000  Decimal  Places”  - 
George  W.  Reitwiesner. 

b.  BRL  TN  381  -  "A  Table  of  Factorial  Numbers  and  Their 

Reciprocals  From  1!  Through  1000 I  to  20  Signifi¬ 
cant  Digits”  -  George  W.  Reitwiesner  -  1951. 

6.  Examination  of  Iterated  Simpson  Numerical  Integration 

Zero  and  first  order  Bessel  Functions  were  obtained 
by  integrating  the  initial  value  problem  y^  =  -  y2 

y2  ■  y^  -  y^x  with  initial  conditions  y^(0)  «  1, 

y2(0)  =  0  using  a  numerical  method  suggested  by 
R.  F.  Clippinger  and  B.  Dimsdale.  The  integration  method 
is  fourth  order  and  has  the  advantage  over  other  methods 
in  that  30-60$  higher  order  systems  can  be  handled  with 
the  same  machine  memory.  Like  the  Runge-Kutta  method, 
it  proceeds  from  initial  data  at  one  point  without  need 
for  data  at  other  points.  Thus  the  grid  size  can  be 
easily  changed  and  machine  storage  requirements  are  not 
unduly  severe.  The  method  is  described  in  Sections  8  and 
8’  of  the  ’’Coding  Notes”  by  Clippinger  and  Dimsdale  and  a 
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ERL  Technical  Note  entitled  "The  Iterated  Simpson  Method 
of  Numerical  Integration"  -  W.  Barkley  Fritz  -  19^2. 

7.  Numerical  Solution  of  Parabolic  Partial  Differential 
Equations 

A  numerical  solution  of  a  linear  parabolic  partial 
differential  equation  with  constant  coefficients  of  the 
type: 

u  +  2  a  u  +  bu  =  k  u. 
xx  x  t 

was  obtained o  A  three  line  difference  equation  approxi¬ 
mation  was  used* 

BRL  769  -  "On  the  Numerical  Solution  of  Parabolic 
Partial  Differential  Equations"  - 
Werner  Leutert  -  195>1° 

8*  Matrix  Inversion 

Inverses  of  matrices  of  orders  £  through  9  using 
partitioning  have  been  obtained  on  the  Eniac*  More 
recently  investigations  on  the  Edvac  and  Ordvac  have 
been  successful  in  inverting  matrices  of  higher  order* 

9*  Maddida  Analogue 

Sin  x  and  cos  x  were  generated  by  a  method  analogous 
to  that  used  by  Maddida  in  order  to  determine  the  ac¬ 
curacy  to  be  expected  by  that  computer*  The  equations 
solved  were  y1  “  z  and  z1  =  -y  using  0*02,  0.002,  and 
0.0002  as  intervals  for  the  independent  variable.  Two 
Eniac  accumulators  were  programmed  similar  to  the  "X" 
and  HR"  registers  of  the  Maddida.  Results  showed  that, 
using  the  smallest  interval  size,  four  significant 
figures  could  be  obtained.  The  study  was  suggested  by 
J.  von  Neumann. 

10*  Characteristic  Trend  Analysis 

The  special  matrix 


was  diagonalized  for  r  »  2  and  n  »  1(1)  3k  as  part  of  a 
problem  in  characteristic  trend  analysis.  More  recently 
the  Edvac  has  solved  the  same  problem  for  r  **  3(1)  -10 
for  the  same  values  of  n. 
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BRL  8U4  -  ”The  Characteristic  Hoots  of  Certain  Matrices”  - 
B.  Dimsdale  and  I.  R.  Hershner  -  1952. 

The  total  machine  time  spent  on  this  group  of  10  problems 
was  of  the  order  of  1*00  machine  hours* 

11.  Miscellaneous  Topics 

Other  mathematical  techniques  were  studied  for  the 
Eniac,  but  were  not  directly  connected  with  any  particular 
machine  program. 

a.  BRL  613  -  ”A  Study  of  Fourth  Order  Interpolation  on 

the  Eniac”  -  H.  B.  Curry  and  M.  Lotkin  -  I9I4.6. 

b.  BRL  632  -  “Inversion  on  the  Eniac  Using  Oscillatory  In¬ 

terpolation”  -  M.  Lotkin  -  19k7  * 

c.  BRL  6l5  -  ”A  Study  of  Inverse  Interpolation  on  the 

Eniac”  -  H.  B.  Curiy  and  W.  'Wyatt  -  19k6, 

II.  Exterior  Ballistic  Laboratory 

Problems  for  the  Exterior  Ballistic  Laboratory  have  been  of  two 
general  classes:  A.  -  Reductions  of  data  from  various  ranges,  and 
B.  -  Theoretical  investigations  in  supersonic  flow  and  related  problems. 

A.  Reduction  of  Range  Data 

1.  Interferometric  Reductions 

This  reduction  program  determines  density  ratios  of 
the  air  surrounding  free  flight  missiles  by  the  reduction 
of  measurements  from  interferograms  obtaining  solutions 
of  the  integral  equation  relating  air  density  within  the 
disturbance  to  the  fringe  shift  on  the  interferogram. 

BRL  797  -  “Interferometric  Analysis  of  Airflow  about  • 
Projectiles  in  Free  Flight”  -  F.  D.  Bennett, 

W.  Co  Carter  and  V.  E.  Bergdolt  -  1952. 

2.  Taw  and  Swerve  Reductions 

The  yaw  reduction  determines  values  of  ten  parameters 
of  the  yawing  motion  of  a  free  flight  missile  by  means  of 
an  iterative  differential  corrections  least  square  pro¬ 
cedure*  The  swerve  reduction  determines  eight  parameters 
by  a  direct  least  square  technique.  AH  data  that  has 
been  reduced  was  obtained  from  photographs  taken  of  missiles 
in  flight  in  the  Aerodynamics  and  Transonic  Firing  Ranges  at  BRL. 


This  problem  has  been  more  successfully  handled  by  the 
Bell  Computers* 

M.I.T.  Department  of  Electrical  Engineering,  Center 
of  Analysis,  Technical  Report  No*  k  - 
“A  Manual  of  Reduction  of  Spinner  Rocket 
Shadowgrams”  -  Z.  Kopal,  K.  E.  Kavanagh, 
and  N.  K.  Rodier  -  19k9. 

B.  Theoretical  Investigations  in  Supersonic  Flow  and  Related 
Problems 

1.  Supersonic  Airflow 

Several  variations  of  the  problem  of  a  theoretical 
treatment  of  supersonic  airflow  about  bodies  of  revolu¬ 
tion  with  and  without  rotation  have  been  completed* 
These  problems  deal  primarily  with  obtaining  solutions 
of  systems  of  hyperbolic  partial  differential  equations. 
A  three  point  Aitken  interpolation  has  been  applied  to 
portions  of  the  Eniac  results  in  order  to  obtain  stream 
data  for  various  Mach  numbers  not  directly  computed  on 
the  machine.  The  following  large  group  of  references 
indicate  the  scope  of  the  work: 

a*  MTAC  III  No.  23  July  19U8  -  “Conputation  of  the 
Airflow  about  a  Cone  Cylinder”  -  M.  Lotkin. 

b.  MTAC  III  No.  23  July  19U8  -  “Airflow  Problem 

Planned  for  the  Eniac”  -  R.  F.  Clipping er. 

c.  BR1  719  -  ” Supersonic  Flow  Over  Bodies  of  Revolu¬ 

tion  (With  Special  Reference  to  High  Speed 
Computing)”  -  R.  F.  Clippinger  and  N.  Gerber  - 

19&. 

d.  BRL  729  and  BRL  730  -  “Tables  of  Supersonic  Flows 

About  Cone  Cylinders*  Part  I:  Surface  Flows; 
Part  II:  Complete  Flows”  -  R.  F.  Clippinger, 

J.  H.  Giese  and  W.  C.  Carter  -  1930. 

e.  BRLM  51k  -  “Theoretical  Supersonic  Pressure  Distri¬ 

butions  on  Non-Yawing  Cone  Cylinders  with  Boat- 
tails”  -  W,  C.  Carter  -  1930. 

f.  BRL  796  -  “Supersonic  Flow  with  Vorticity  About  a 

Slightly  Yawing  Body  of  Revolution”  - 
W.  C.  Carter  -  1932. 
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g.  BRL  813  -  "On  the  Numerical  Solution  of  Hyperbolic 
Systems  of  Partial  Differential  Equations  with 
Two  Characteristic  Directions"  -  W.  C.  Carter 
and  G.  I.  Spencer  -  1952. 

2.  Supersonic  Taylor-Maccoll  Flow  Past  Non-Yawing  Cones 

A  system  of  three  ordinary  differential  equations  was 
solved  to  determine  the  Taylor-Maccoll  flow  past  a  non¬ 
yawing  cone.  The  solution  was  obtained  as  a  characteristic 
curve  with  the  components  of  the  velocity  of  air  along  the 
curve.  The  density  and  pressure  of  the  air  were  determined 
from  the  velocity  and  terminal  conditions  and  these  density 
distributions  and  interferometer  fringe  shifts  were  cal¬ 
culated  with  y  *  for  cones  with  semi-vertex  angles  of 

5  K°,  K  “  2(1)8,  over  the  Mach  number  range  of  1.3,  1.7* 

2.0,  3.0#  li.O,  and  5*0. 

BRL  793  -  "Densities  and  Interferometric  Fringe  Shifts 
for  Taylor-Maccoll  Flows"  -  J.H.  Giese  -  193>2. 

3.  Prandtl-Meyer  Flow 

A  table  of  values  of  density,  temperature,  and  pres¬ 
sure  ratios  of  supersonic  Prandtl-Meyer  Flow  was  confuted. 
The  equation: 

<La!  . 

d  e 

was  integrated  by  the 

6  =  0°  (0.05)  102.50°. 

k*  Mathieufs  Equation 

Solutions  were  obtained  of  the  Mathieu  Equation: 

y"  +  1/1;  (h  +  ycos  x)  y  *  0  for  two  sets  of  initial 
conditions : 

a)  7  -  1,  7*  -  0 

t>)  y  88  °>  y*  s  1 

The  method  of  Runge-Kutta  was  used  to  integrate  the 
equation  from  x  *  0  to  x  -  n  listing  values  at  0.1  radian 
intervals. 

The  range  of  parameters  follow: 

2 

n  :  0.2  (0.2)  12. 0;  the  interval  in  y  varied 

from  0.25  for  small  y  to  2.0  for  large  y. 


im 


method  of  Iterated  Simpson  for 
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Listings  of  the  results  are  available  from  the  Com¬ 
puting  Laboratory,  BRL. 

Earlier  Mathieu  Equation  solutions  were  obtained  by 
the  Moore  School  of  Electrical  Engineering  of  the  University 
of  Pennsylvania. 

5*  Two-dimensional  Flexible  Throat  Wind  Tunnel  Nozzles 

The  problem  had  as  its  object  the  determination  of  two- 
dimensional  nozzle  shapes  for  a  flexible  throat  wind  tunnel. 
The  results  were  used  to  determine  cam  shapes  for  an  auto¬ 
matic  jack  control  system  for  test  section  Mach  numbers 
ranging  from  1.25  (0.25)  5*50.  The  problem  involves 
selection  of  suitable  combinations  of  three  parameters 
to  produce  nozzles  which  have  the  prescribed  test  section 
Mach  numbers  and  which  satisfy  three  additional  conditions 
imposed  by  the  physical  characteristics  of  the  tunnel  and 
by  the  desire  to  have  reasonable  cam  shapes.  Since  the 
choice  could  not  be  determined  readily  by  theoretical 
analysis,  a  trial  and  error  procedure,  requiring  several 
computations  for  each  nozzle,  was  employed. 

6.  Shock  Curvature  vs  Nose  Curvature 

An  extension  of  some  work  by  C.  C.  Lin  relating  to  a 
comparison  of  shock  curvature  and  body  curvature  about 
the  nose  of  pointed  objects  in  supersonic  flow. 

III.  Interior  Ballistic  Laboratory* 

A.  Interior  Ballistic  Trajectories 

Various  Interior  Ballistic  Tables  for  powder  of  constant  burning 
surface  based  on  an  interior  ballistic  system  described  by  Laidler  and 
Hershey  have  been  computed.  The  integration  methods  of  Kutta-Simpson 
and  Runge-Kutta  have  been  utilized  for  various  areas  of  the  tables.  Ex¬ 
tensive  additional  computations  have  recently  been  completed  on  the  Ordvac. 

BRIM  552  -  "Interior  Ballistic  Integrations  for  Powder  of 
Constant  Burning  Surface"  -  S.  Kravitz  -  19S>1  • 

B.  Mathematical  Theory  of  Thermal  Ignition 

A  theoretical  investigation  of  ignition  and  burning  of  powder 
by  obtaining  solutions  of  a  non-linear  parabolic  partial  differential 
equation  using  six  grid  sizes  to  aid  in  analyzing  error. 

.  — l/u 

Uj  -  U££  +  e  ' 

in  the  semi-infinite  region  £  >  0  . 

*  See  also  VI,  B,  1 
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The  boundary  conditions  are 

u  ^  =  -H  (u  g  -  ^  £=0,  T  -  0 

u  -*U  as  <$  — *  *• 

Initial  conditions  are  u  =  U,  £  >  0,  and  T  *  0, 

BRL  756  -  ’’Mathematical  Theory  of  the  Ignition  Process 
Considered  as  a  Thermal  Reaction"  - 
B.  L.  Hicks,  J.  W.  Kelso  and  J.  Davis  -  1951* 

C.  Crystal  Structure 

A  numerical  method  for  the  determination  of  the  Crystal 
Structure  of  Banfield's  or  Kenyon's  free  radical  was  investigated. 

The  numerical  problem  was  to  determine  a  sum  of  22  products  of  an  ex¬ 
pression  involving  sines  and  cosines  of  the  form  cos  2n  H  x^  • 

cos  2ti  K  y^  *  sin  2tt  L  z^  where  H,  K,  and  L  are  positive  integers 

and  y\ ,  z^  are  the  space  coordinates  of  the  22  carbon,  nitrogen 

and  oxygen  atoms  in  Kenyon's  radical  obtained  by  the  methods  of 
crystallography. 

There  were  U  sums  for  each  set  of  the  integers,  the  exact  form 
of  the  product  being  determined  by  whether  K  is  odd  or  even.  In  all, 
819  cases  involving  the  computation  of  over  100,000  sines  or  cosines 
were  completed.  The  problem  was  also  used  as  a  trial  problem  to  gain 
experience  with  the  operation  of  the  Uni  vac  of  the  Eckert-Mauchly  Com¬ 
puter  Corporation. 

IV.  Ballistic  Measurements  Laboratory 

The  main  problem  handled  for  BML  has  been  the  reduction  of  Doppler 
Data  from  three  or  five  simultaneous  radio-Doppler  observations.  The 
program  has  been  used  to  reduce  data  for  about  20  firings  to  date  and 
determines r 

A.  rectangular  coordinates  and  probable  error, 

B.  velocity  and  acceleration  by  numerical  differentiation,  and 

C.  reduction  of  rectangular  coordinates  to  transverse  mercator 
coordinates  of  the  missile. 

The  method  of  least  squares  is  used.  The  three  station  problem 
is  discussed  in  the  first  reference,  the  five  station  problem  in  the 
second,  and  additional  information  is  contained  in  the  third. 
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1*  BEL  638  -  "Doppler  Determination  of  Position"  -  Boris  Garfinkel 

19h7o 

2.  BEL  707  -  "Least  Square  Determination  of  Position  from  Radio 
Doppler  Data"  -  Boris  Garfinkel  -  19li9« 

3«  MTAC  HI  No.  25  January  1 9k9  -  "A  Comparison  of  Various  Com¬ 
puting  Machines  Used  in  the  Reduction  of  Doppler 
Observations"  -  Dorrit  K off lei t. 


V.  Surveillance  Laboratory 

Although  several  problems  were  originally  initiated  by  this  group 
including  the  Characteristic  Trend  Analysis  Problem,  the  only  one  we 
shall  discuss  in  this  section  is  one  in  probability  distribution.  This 
problem  had  to  do  with  the  probability  distribution  of  the  difference 
between  the  largest  observation  and  the  sample  mean  (or  the  difference 
between  the  sample  mean  and  the  smallest  observation)  for  normally  dis¬ 
tributed  random  variates.  The  table  of  probabilities  for  these  sta¬ 
tistics  is  very  useful  in  determining  whether  the  greatest  or  smallest 
observations  in  a  sample  of  n  observations  can  be  judged  "outlying" 
observations.  For  normally  distributed  variables  which  are  ranked  after 
the  sample  is  drawn,  i.  e.  into  the  form  u^  6  Ug  £  . ..  £  ufl,  then  the 

functions  Hn  (nk)  give  the  probability  that  the  extreme  deviation  from 

the  sample  mean,  i.e.  uQ  -  u  or  u  -  u^,  will  not  exceed  k  times  the 

standard  deviation  of  the  normal  population  from  which  the  random  sample 
of  size  n  was  drawn.  The  functions  Hn  (k)  are  defined  below  and  were 

computed  on  the  Eniac  for  n  =  2(1)  2$ 


fit  -  te)  <*w 

o 

Tables  of  the  above  functions  which  have  use  not  only  for  the 
problem  of  outlying  observations  but  also  for  the  more  general  field 
of  Order  Statistics  have  been  published  in  the  March  1950  issue  of  the 
Annals  of  Mathematical  Statistics,  pp  27-58.  "Sample  Criteria  for 
Testing  Outlying  Observations"  -  Frank  E.  Grubbs. 


H1  >) 


=  1 


Hn  ^  "  /  n 
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VI.  Outside  Groups 

A.  Atomic  Energy  Commission 

Eight  problems  have  been  completed  for  the  A.  E.  C.  The  informa¬ 
tion  given  is  limited  to  unclassified  reports  published  by  the  individual 
sponsors  of  the  particular  problems. 

1.  Calculations  in  the  Liquid-Drop  Model  of  Fission 

The'  liquid  drop  model  of  fission  as  developed  by  Bohr, 
Wheeler,  and  others  was  investigated.  Classical  deforma¬ 
tion  energies  were  calculated  for  many  shapes  not  all  near  * 
spherical.  The  shapes  considered  were  axially  and  primarily, 
but  not  exclusively,  bilaterally  symmetric,  and  single 
valued  in  polar  representations.  Agreement  was  obtained 
with  observed  fission  thresholds  and  spontaneous  fission 
rates  obtained* 

"Calculations  in  the  Liquid-Drop  Model  of  Fission"  - 
S.  Frankel  and  N.  Metropolis  -  Physical 
Review  Vol.  72  No.  10,  Nov.  19l;7>  PP  911+-925* 

2.  Fermi -Thoraas-Dirac  Equation 

Solutions  of  the  Fermi-Thoraas-Dirac  equation: 

-iX-  «  x  (  e  ♦  yVZ  *'1/2)  3 

d  "2L 

-3  _c!  _p  _2 

where  «  3  •  2  -n  z 

were  obtained  in  terms  of  the  variable  w  **  (2x)"^  ♦  2  yr 
is  tabulated  to  5  decimal  places  at  intervals  of  0.08 
extending  until  ~[j/  becomes  negative.  The,  parameter  z 
takes  on  2k  values  6(U)  lit,  16,  18(U)  26,  29(lt)  81,  - 
8U(U)  92  and  there  are  8  different  initial  slopes  . 

"Solutions  of  the  Fermi-Thomas-Dirac  Equation"  - 

N.  Metropolis  and  J.  R.  Reitz,1  Journal  of 
Chemical  Physics  Vol.  19,  1951,  PP*  555-573. 

3*  Table  of  Atomic  Masses  -  Program  ran  on  card  control. 

"Table  of  Atomic  Masses"  published  by  Argonne 
National  Laboratory  -  March  1950. 
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U.  Miscellaneous  Problems 

The  following  references  indicate  the  scope  of  further 
unclassified  work: 

a*  Seminar  on  Scientific  Computation  -  Nov*  19h9  -  IB1!  Corp. 
"Modification  of  the  Monte  Carlo  Technique"  - 
H0  Kahn* 

b*  Proceedings  of  the  Monte  Carlo  Conference  -  "Estimation 
of  Particle  Transmission.  Through  Plane  Slabs"  - 
T.  E.  Harris  and  H*  Kahn. 

c.  Journal  of  Applied  Physics  Vol.  21  No.  3,  PP  232-237  - 
March  195>0  -  "A  Method  for  the  Numerical  Calcula¬ 
tion  of  Hydrodynamic  Shocks"  -  J.  von  Neumann  and 
R*  Do  Richtmyer. 

d*  Nat.  Bur«  Stand.  Appl.  Math.  Ser.  No,  12,  36-38  (19i?l)  - 
"Various  Techniques  Used  in  Connection  with  Random 
Digits"  -  J.  von  Neumann.  This  is  a  paper  pre¬ 
sented  at  the  Monte  Carlo  Method  Symposium 
describing  a  method  of  producing  random  numbers 
according  to  an  assigned  probability  law  on  the 
Eniac,  using  random  numbers  from  a  rectangular 
distribution  in  the  interval  (0,1). 

e.  Nat*  Bur.  Stando  Appl.  Math,  Ser.  No.  12,  19-20  (195>1)  - 
"Report  on  a  Monte  Carlo  Calculation  Performed 
with  the  Eniac"  -  M,  Mayer.  Here  the  absorption 
of  neutrons  in  a  composite  cylindrical  body  was 
investigated. 

Bo  U*  S.  Bureau  of  Mines 

1.  Equilibrium  Composition  and  Thermodynamic  Properties 

This  problem  has  as  its  object  the  determination  of  the 
mole  fractions  of  the  gaseous  constituents  of  the  equilibrium 
composition  of  the  four  component  system  of  hydrogen,  oxygen, 
nitrogen,  and  carbon  at  temperatures  of  1000°  K  to  £000°  K 
and  pressures  of  0.1  atmospheres  to  those  of  interior  bal¬ 
listics.  In  addition  various  thermodynamic  properties  have 
been  evaluated.  Approximately  7%  of  the  Eniac* s  time  has 
been  spent  on  this  problem  and  the  work  is  continuing. 

a.  Chemical  and  Engineering  News,  American  Chemical  Society, 
Vol.  27,  5  Sept,  19k9°  "Combustion  Gases  ... 
Equilibrium  Composition  and  Thermodynamic  Proper¬ 
ties"  -  Stuart  R,  Brinkley,  Jr.  and  Bernard  Lewis. 
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b.  Bureau  of  Mines  Report  of  Investigations  I4.806  - 
nThe  Thermodynamics  of  Combustion  Gases: 

General  Considerations"  -  Stuart  R.  Brinkley, Jr. 
and  Bernard  Lewis  -  U.S.  Department  of  the 
Interior  -  April  195>2. 


2.  Optimum  Spacing  and  Performance  of  Gas  Wells 


The  non-linear  parabolic  equation: 


at 


was  solved  for  specific  boundary  conditions  relating  to 
gas  well  performance  and  optimum  spacing  of  such  wells. 

A  single  solution  has  been  obtained  using  a  six  point 
finite  difference  procedure.  Additional  machine  studies 
are  contemplated  on  Edvac  and  Ordvac. 

3.  Well  Depletion 

A  preliminary  study  was  made  of  the  rate  of  withdrawal 
of  oil  from  a  well  consisting  of  two  pressure  areas.  The 
method  of  Runge-Kutta  was  used  to  integrate  three  ordinary 
differential  equations  in  one  independent  variable.  Ad¬ 
ditional  computations  have  been  completed  on  the  Edvac. 

C.  Universities  (including  Government  Contractors) 

1.  Compressible  Laminar  Boundary  Layer  (University  of  Cambridge) 

An  examination  was  conducted  by  D.  R.  Hartree  of  com¬ 
pressible  laminar  boundary  layer  flow  for  the  case  of  a 
flat  plate  at  zero  incidence. 

The  mathematical  method  is  described  in  the  Philo¬ 
sophical'  Transactions  of  the  Royal  Society  of1  London  - 
Series  A  -  Math  and  Physical  Sciences  No.  827,  Vol.  2ljl, 
pp  I-69,  22  June  19 W  -  "The  Laminar  Boundary  Layer  in 
Compressible  Flow"  -  W.  F.  Cope  and  D.  R.  Hartree. 

A  briefer  discussion  is  given  in  MTAC  III  No.  23, 

July  19^8  -  "Laminar  Boundary  Layer  Flow  in  a  Compres¬ 
sible  Fluid"  -  J.  V.  Holberton. 

2*  Zero  -  Pressure  Properties  of  Diatomic  Gases  -  (University 
of  Pennsylvania) 

J.  A«  Goff,  Dean  of  the  Towne  Scientific  School, 
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determined  the  zero-pressure  properties  of  certain  diatomic 
gases  using  the  equations? 

(r“t)  “  +  B//a 

'  7  p  *  o 

-  3*3>  +  x  "  (B/a)2 
p  *  o  a 

where  B  =  T  and  y  *  ^  ^  £ 

d  T  d  T^ 

The  computations  were  based  on  the  best  available 
spectroscopic  data* 

An  additional  problem  for  the  University  of  Pennsylvania 
is  discussed  in  I,  C,  1* 

3,  Refraction  of  Plane  Shock  Waves  (institute  for  Advanced 
Study) 

Numerical  solutions  were  obtained  of  the  equations 
governing  the  position  of  the  various  relevant  quantities 
in  an  assumed  shock  configuration  for  the  refraction  of 
plane  shocks.  The  assumed  shock  configuration  consisted 
of  an  incident  shock s  a  reflected  shock,  and  a  transmitted 
one.  It  was  further  assumed  that  the  pressure  is  constant 
in  the  angular  domains  between  these  shocks  and  that 
across  the  density  discontinuity  responsible  for  the  re¬ 
fraction,  the  pressure  is  continuous  as  is  the  deflection 
of  the  flow. 

a.  “Reflection  of  Plane  Shock  Waves11  -  A*  H.  Taub, 

Physical  Review,  Vol.  72,  pp  f>l-60  (19li7)» 

b.  “Some  Numerical  Results  on  Refraction  of  Plane 

Shocks'*  -  A.  Ho  Taub  -  University  of 
Washington. 

U.  AxLsymmetric  Wind  Tunnel  Nozzle  Design  (M.  I.  T«) 

Determination  of  axisymmetric  wind  tunnel  nozzles 
solving  a  system  of  five  first  order  partial  differential 
equations. 

BRL  79^  -  "Supersonic  Axially  Symmetric  Nozzles"  - 
R.  F.  Clippinger  -  195>2« 


2h 


5.  Mie  Theory  of  Light  Scattering  Functions  for  Spherical 
Particles  (University  of  Michigan) 

Evaluation  of  electromagnetic  wave  scattering  functions 
for  spherical  particles  to  be  used  in  connection  with  ex¬ 
perimental  work  in  the  examination  of  small  particles  ac¬ 
cording  to  the  Mie  Theory.  Two  books  have  appeared  giving 
the  Eniac  results: 

a.  "Tables  of  Riccati  Bessel  Functions  for  Large  Arguments 

and  Order",  and  "Light  Scattering  Functions  for 
Spherical  Particles  Part  I,  Part  II", •  both  by 
R.  0*  Gumprecht  and  C.  M.  Sliepcevich,  Engineering 
Research  Institute,  University  of  Michigan,  1951* 

b.  Journal  of  the  Optical  Society  of  America,  Vol  h2. 

No.  I;,  pp  226-231,  "Angular  Distribution  of 
Intensity  of  Light  Scattered  by  Large  Droplets 
of  Water"  -  R.  0*  Gumprecht,  Neng-Lun  Sung, 

Jin  H.  Chin,  and  C.  M.  Sliepcevich  -  University 
of  Michigan  -  April  1952. 

6.  Investigations  in  Weather  Forecasting  (Institute  for  Ad¬ 
vanced  Stucty) 

A  method  was  adapted  to  the  Eniac  for  the  numerical 
solution  of  the  barotropic  vorticity  equation  over  a 
limited  area  of  the  earth* s  surface. 

"Numerical  Integration  of  the  Barotropic  Vorticity 
Equation"  -  J.  G.  Charney,  R.  Fj or toft, 

J.  von  Neumann  -  Tellus,  Vol.  2,  No.  1;  - 
November  1950* 

The  problem  discussed  in  I,  C,  3  was  performed  by 
D.  H.  Lehmer  of  the  University  of  California. 

D.  Other  Outside  Agencies 

1.  Heat  Flow  (Armour  Research  Institute) 

A  group  of  96  one  dimensional  heat  flow  problems  with 
a  single  source  and  a  new  set  of  boundary  conditions  for  each 
problem  was  computed.  .  An  error  study  and  16  of  the 
one  and  two  slab  cases  were  successfully  completed  before 
the  project  was  cancelled. 

2.  Guided  Missile  Trajectories  (General  Electric  Co.) 

Two  problems  in  guided  missile  trajectory  computations 
have  been  completed  on  the  Eniac. 
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a.  Power  flight  computation  for  the  two  dimensional 
case  with  cartesian  coordinates,  a  time  and  posi-  1 
tion  dependent  control  system,  the  moment  equation, 
and  a  time  and  atmosphere  pressure  dependent 
thrust.  The  integration  of  all  quantities,  es¬ 
sentially  a  sixth  order  differential  equation, 

was  computed  by  the  iteration  of  the  trapezoidal 
rule. 

b.  Free  flight  trajectories  for  the  two  dimensional 
case  with  intrinsic  coordinates,  i.  e*  velocity 
and  path  angle,  and  a  control  system  determined 
by  a  second  order  differential  equation  and  by  a 
time  dependent  signal.  The  integration  of  the 
three  second  order  derivatives  was  computed  by 

an  iterated  use  of  the  trapezoidal  rule  and  first 
derivatives  by  the  Euler  MacLaurin  rule. 


W. 


Barkley  Fritz 
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